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DYNAMICAL  SYSTEMS  IN  NEUROPHARMACOLOGY 

FINAL  REPORT 

MARK  A.  GEYER,  Ph.D. 

AND 

ARNOLD  J.  MANDELL,  M.D. 

The  original  proposal  for  this  project  listed  three  principal  goals: 

1.  Further  elaboration  and  application  of  the  quantitative  techniques  of  modem  dynamical  systems 

theory  and  experimental  mathematics  to  sequential  data  or  "time  series". 

2.  Numerical  exploration  of  both  smooth  and  discrete  abstract  dissipative  dynamical  systems  as  real¬ 
istic  low-dimensional  compressions  of  the  dynamical  response  pattern  to  increasing  doses  of 
stimulants. 

3.  Development  of  phase-space  gnq)hic  computer  programs  for  the  qualitative  desciipdcm  of  changes 
in  complexity  induced  by  stimulant  and  other  drugs  in  r  oronal,  neurochemical,  and  behavioiai 
dynamical  systems. 

Data  were  to  be  generated  in  three  laboratories  studying  the  effects  in  rats  of  stimulants  and 
other  parametirci  perturbations  on:  (a)  locus  coetuleus,  substantia  nigra,  and  cortical  neuron  interspUce 
intervals;  (b)  continuously  HPLC(ec)-monitored  levels  of  NE,  MHPG,  DA,  DOPAQ  HVA,  5HT,  and 
SHIAA  in  brainstem  perfusion  studies  of  freely  moving  rats;  (c)  digitally  resolved  temporal  and  ^- 
tial  exploratory  behavior  in  conjunction  with  (b).  The  laboratories  to  be  involved  were  diose  of  Drs. 
Stephen  Foote  (a),  Ronald  Kuczenski  (b),  and  David  Segal  (c)  of  UCSD.  The  shared  interests  of 
these  investigators  in  brainstem  biogenic  amine-mediated  neurobiological  mechanisms  of  bduivior 
and  the  use  of  cocaine,  amphetamine,  and  other  pharmacological  probes  of  mechanisim  was  expected 
to  permit  cross-level  exploration  for  the  dynamical  univosality  pr^cted  by  the  theory. 

As  the  project  progressed,  some  alterations  in  the  emphasis  of  the  project  were  required.  The 
fundamental  premise  of  the  project  was  to  assess  the  utility  of  applying  non-linear  dynamicai 
methods  to  high-quality  biological  data.  That'  general  objective  of  the  project  was  itecoiiq)lished. 
However,  the  specific  databases  upon  which  these  method  were  proposed  to  be  develtqied  were  not 
all  used  effectively.  Specifically,  the  planned  use  of  time  series  data  derived  fixxn  Dr.  Kuczensld’s 
HPLC  analyses  of  neurochemical  changes  in  freely  moving  rats  never  materialized.  Similariy,  ik> 
substantive  progress  was  made  using  data  from  Dr.  Segal's  locomotor  activity  chambers.  Sow 
promising  work  has  been  accomplished  in  collaboration  with  Dr.  Foote,  using  his  measures  of  single 
unit  responses  of  cells  in  the  locus  coeruleus,  substantia  nigra,  and  cortex  of  rats.  This  wrak  is 
currendy  in  manuscript  form  and  should  be  submitted  for  publication  within  the  year. 

In  some  of  our  more  stricdy  theoretical  wtnk,  we  demonstrated  some  dynamical  consequences 
of  making  the  fomnal  characteristic  of  neuronal  elements  and  their  connectivities  nxxe  realistic  than 
those  of  modem  neural  computer  algorithms  [3].  Specifically,  we  studied  two  multiplicatively  cou¬ 
pled  logistic  maps,  each  with  the  quadratic  response  characteristics  of  neurons  widi  autoiecqMor 
mechanisms,  which  iteratively  exchange  their  outputs.  The  dynamic  phenomenology  of  this  two- 
netnon  system,  though  noninvertible,  demonstrated  characteristics  of  both  automoiphisms  of  die  real 
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line  and  diffeoiiKnphisnis  of  the  plane,  with  snHX>th  parametric  transitions  between  them.  We  por¬ 
trayed  a  unified  parameter  space  containing  the  negative,  positive,  and  mixed-coupled  regimes, 
categorized  the  phase  space  structures  using  characteristic  exponent  inequalities,  expkned  rqnesrata- 
tive  parametric  paths,  demonstrated  generic  bifurcation  sequences,  and  reported  number  thetmtic  ord¬ 
ering  in  parameter  space.  We  conjectured  that  most,  if  not  all,  "middle-layer-like"  brain  systems,  in 
contrast  to  primary  sensory  and  motor  information  transport  systems,  are  dominated  by  intrinsic 
dynamics  of  the  sort  we  demonstrated.  According  to  this  view,  external  input  serves  as  a  perturba¬ 
tion  of  these  already  ongoing  complex  systems.  The  intrinsic  instabilities  of  these  middle-layer 
dynamics  return  the  unpredictability  required  by  theroetical  studies  of  real  brain  functions  to  noodels 
of  neural  networks.  We  demonstrated  a  variety  of  parameter-dependent  intrinsically  ordered  attractor 
sequences  and  number  theoretic  series  which  emerge  fiom  the  global  dynamics  and  exemplify  the 
rich  dynamical  "machine"  language  potential  for  specific  coding  by  brain-like  information  transport 
and  learning  algorithms. 

Another  theoretical  development  was  pursued  to  examine  the  statistical  properties  of  real  macro- 
molecular  sequences  using  an  approximation  based  on  the  composition  of  two  ergodic  sets  [2].  Both 
varied  and  reliable  mechanisms  of  information  transport  in  molecular  biological  systems  can  be 
obtained  fiom  the  universality  properties  of  continuous  and  discrete  dynamical  systems  whose  output 
has  the  general  properties  of  a  set  which  combines  sensitivity  to  initial  conditions  and  indeconq)osa- 
bility  with  a  positive  density  of  periodic  points.  This  approach  required  a  revision  of  infcmnatitm 
theory  for  macromolecular  systems  dictated  by  their  failure  to  fulfill  the  conditions  of  Shaimon’s 
existence  proof:  there  exists  a  one-to-one  mapping  that  makes  the  channel  capacity  finite.  Relaxa¬ 
tion  studies  of  polypeptides,  proteins,  and  nucleotides  demonstrate  a  lack  of  characteristic  times. 

Using  calcitonin  as  the  exemplar,  a  more  appropriate  statistical  approach  to  macromolecular  coding 
was  explored  using  a  hyperbolic  distribution  with  an  infinite  multiplicity  of  moments,  multiplicative 
noise,  and  belonging  to  the  family  of  convolutionally  stable  distributions  without  finite  mean  vari¬ 
ance  [2].  This  work  used  a  phase  space  reconstruction  technique  to  relate  the  potencies  of  impratant 
calcitonin  analogues  to  the  global  "tightness"  of  the  vorticity. 

The  most  significant  work  related  to  concrete  biological  systems  was  based  on  die  data  from  Dr. 
Geyer’s  Behavioral  Pattern  Monitor  (BPM)  system,  which  provided  the  pilot  data  on  behavioral 
measures  presented  in  the  original  proposal.  The  extensive  database  derived  from  the  BPM  system 
proved  to  be  ideal  for  addressing  the  conceptual  questions  originally  proposed.  Thus,  we  had  consid¬ 
erable  success  by  carefully  examining  data  from  studies  of  the  behavioral  profiles  of  psychoactive 
drugs  as  characterized  in  the  BPM  system.  In  particular,  we  developed  new  methods  based  on  pro¬ 
cedures  developed  in  the  context  of  theoretical  physics  and  non-linear  dynamics.  Using  a  k-D  split¬ 
ting  algorithm,  we  described  the  initial  demonstrations  of  the  utility  of  complexity  and  entropy  meas¬ 
ures  in  biological  applications,  using  data  from  the  BPM  system  [7].  Measures  of  complexity  derived 
from  ergodic  theory  of  dynamical  systems  were  developed  and  applied  to  an  exemplary  data  set 
describing  locomotor  movements  of  rats  in  a  bounded  space.  A  symbolic  dynamical  system  was 
obtained  by  partitioning  the  event  space  into  equally  probable  partition  elements,  using  a  k- 
dimensional  tree.  The  measures  calculated  from  the  symbolic  sequences  included  the  topological 
entropy  (ht)-i.e.,  the  rate  of  increase  of  all  possible  sequences  with  increasing  sequence  length-and 
the  metric  entropy  (hm)-i.e.,  the  rate  of  increase  of  all  likely  sequences  with  increasing  sequence 
length.  These  measures  were  used  to  assess  changes  in  rat  locomotor  behavior  as  recorded  in  the 
behavioral  pattern  monitor  (BPM)  that  were  induced  by  amphetamine  (0.25,  0.50,  1.0,  (x  2.0  mg/kg) 
and  3,4-methylenedioxymethamphetamine  (MDMA;  1.25,  2.5,  5.0,  or  10.0  mg/kg).  Amphetamine 
increased  the  mean  activity,  ht,  and  hm.  MDMA  resulted  in  a  monotonic  dose-re^xmse  curve  fiv 
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activity  but  exhibited  a  biphasic  dose  response  in  ht  and  hm.  In  particular,  some  animals  in  the 
higher  dose  groups  showed  a  ht  in  the  range  of  the  saline  controls,  whereas  other  animals  exhibited  a 
significantly  reduced  ht  and  a  greater  decrease  in  hm,  suggesting  that  two  different  behavioral  reac¬ 
tions  coexist  within  the  same  higher  dose  range  of  MDMA.  An  important  implication  of  our  method 
is  that,  in  applied  ergodic  measure-theoretic  approaches,  the  partition  that  determines  the  elements  of 
the  symbolic  dynamical  system  should  not  be  specified  a  priori  on  abstract  mathematical  grounds  but 
should  be  chosen  relative  to  its  significance  widi  respect  to  the  data  set  in  question.  Here,  the  animal 
constructs  its  own  spatiotemporal  partition  in  behavioral  phase  space. 

In  addition  to  the  complexity  and  entropy  measures,  we  also  made  great  progress  in  the  use  of 
scaling  measures  to  provide  more  sensitive  characterizations  of  drug  effects.  We  used  scaling  meas¬ 
ures  to  provide  more  sensitive  characterizations  of  drug  effects  [4,6].  Conceptually,  these  measures 
are  derived  from  the  concept  of  dimensionality  in  the  field  of  firactal  geometry.  We  have  defined  two 
global  scaling  measures  fix>m  our  BPM  data,  one  characterizing  the  temporal  structure  of  the  animal’s 
behavior,  the  other  being  sensitive  to  changes  in  the  geometrical  structure  of  the  animal’s  locomotcnr 
path.  These  two  measures  enable  us  to  assess  quantitatively  similarities  and  differences  between 
different  drug  treatments.  The  first  full  description  of  this  approach  appeared  in  Psychopharmacology 
[6].  This  paper  details  the  derivations  and  applications  of  the  temporal  scaling  exponent,  alpha,  and 
the  spatial  scaling  exponent,  d,  and  introduces  the  d-alpha  plane.  Basically,  alpha  assesses  the  level 
of  activity  of  an  animal  and  has  proven  to  be  more  consistent  and  sensitive  than  the  traditional  meas¬ 
ure  of  photobeam  breaks,  presumably  because  it  utilizes  all  the  available  information.  The  spatial 
scaling  exponent  d  describes  the  roughness  of  the  path  traversed  by  the  animal  -  it  assesses  the  aver¬ 
age  geometrical  information  contained  in  locomotor  activity  in  a  maimer  that  is  independent  of 
changes  in  the  amount  of  activity.  The  d-alpha  plane  provides  a  condensed  way  to  present  and  visu¬ 
alize  a  large  amount  of  information  about  the  nature  of  a  drug’s  effects  across  a  range  of  doses  and 
facilitates  comparisons  between  the  effects  of  different  drugs.  To  illustrate  this  feature  of  our 
approach,  we  compared  the  effects  of  a  variety  of  psychostimulant  drugs  and  used  these  comparisons 
to  address  hypotheses  regarding  the  behavioral  processes  affected  by  these  drugs  [4,5,6]. 

A  variety  of  psychoactive  substances  (amphetamine,  nicotine,  scopolamine,  apomorphine, 
lisuride,  and  MDMA)  were  tested  to  examine  whether  a  scaling  hypothesis  is  appropriate  for  the 
description  of  the  amount  and  the  structure  of  rat  locomotor  paths  recorded  in  the  BPM  [6].  The 
analytical  approach  was  based  on  the  assumption  that  the  scaling  behavior  of  a  few  collective  vari¬ 
ables  may  characterize  sufficiently  changes  in  the  animal’s  behavior  induced  by  different  drugs.  The 
temporal  scaling  exponent  alpha,  describing  the  ratio  of  fast  to  slow  responses  in  the  BPM,  sensi¬ 
tively  reflected  the  different  stimulant  properties  of  the  substances.  The  spatial  scaling  exponent  d, 
which  relates  the  average  path  length  to  the  resolution  used  to  irteasure  consecutive  responses,  was 
found  to  discriminate  substances  that  had  been  separated  previously  via  qualitative  descriptions. 
Several  behavioral  response  categories  emerged  from  comparisons  of  the  locations  of  different  drugs 
on  a  two-dimensional  d-alpha  plane.  Scopolamine,  MDMA,  lisuride,  and  high  doses  of  apomorphine 
increased  alpha  while  decreasing  d,  whereas  amphetamine,  nicotine,  and  caffeine  produced  an 
increased  alpha  with  no  change  or  an  increase  in  d.  Stereotypies  could  be  identified  on  the  opposite 
ends  of  the  spatial  scaling  exponent  scale  and  were  interpreted  as  reflecting  two  kinds  of  persevera¬ 
tion.  These  results  suggest  that  scaling  approaches  can  used  to  assess  quantitatively  die  state  of 
the  animal  based  on  its  locomotor  behavior  and  that  the  exponents  can  serve  as  collective  variables 
providing  a  macroscopic  description  based  on  the  microscopic  elements  of  behavior. 

Using  both  our  traditional  measures  and  the  newer  theoredcally  driven  measures,  we  explored 
the  mechanisms  contributing  to  the  locomotor  hyperactivity  induced  by  MDMA  [1,5].  We 
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discovered  that  the  profound  locomotor  activation  producni  by  MDMA  and  related  drugs  is  ihh  due 
primarily  to  dopamine  release,  but  is  due  to  serotonin  release  [1].  We  blocked  MDMA-induced 
hyperactivity  with  the  selective  serotonin  uptake  blockers  fluoxetine,  sertraline,  or  zimelidine  or  with 
the  serotonin  synthesis  inhibitor  PCPA.  By  contrast,  these  pretreatments  either  potentiated  or  had  no 
effect  on  the  hyperactivity  induced  by  the  dopamine  releasing  agent  amphetamine.  Conversely,  the 
catecholamine  synthesis  inhibitor  AMPT,  which  blocks  the  effects  of  amphetamine,  did  not  a^ect  the 
response  to  MDMA.  We  applied  our  scaling  and  complexity  measures  to  this  problem  [5].  The  use 
of  scaling  measures  enabled  us  to  further  differentiate  the  effects  of  MDMA  from  those  of  either 
amphetamine  or  hallucinogens.  In  addition,  comparisons  of  the  behavioral  effects  of  MDMA  across 
both  repeated  tests  and  repeated  administrations  of  the  drug  confirmed  that  our  spatial  scaling 
exponent,  d,  was  affected  differently  from  more  traditional  measures.  Furthermore,  d  proved  to  be 
the  most  sensitive  of  our  various  measures  of  the  effects  of  serotonin  releasers  on  the  spatial  structure 
of  the  locomotor  paths  exhibited  by  the  animals.  The  combination  of  our  new  and  old  measures  of 
drug  effects  on  behavior  clearly  provides  more  information  than  can  be  obtained  by  reliance  on  any 
one  set  of  measures. 

The  d-alpha  plane  was  also  used  to  systematically  compare  various  methylenedioxy-substituted 
phenethylamines  having  varying  potencies  as  releasers  of  dopamine  or  serotonin  [5].  This  analytical 
approach  enabled  us  to  make  discriminations  between  drugs  and  specific  enantiomers  of  drugs  that 
have  previously  required  much  more  elaborate  and  time-consuming  procedures  such  as  drug- 
discrimination  tests.  With  this  approach,  we  can  characterize  and  quantify  different  drug  effects  in  a 
single  test  session  using  a  procedure  that  does  not  require  repeated  administrations  of  the  drug  to  the 
animals.  Thus,  the  effects  of  chronic  versus  acute  exposures  to  the  drugs  can  now  be  assessed  reli¬ 
ably. 

In  further  work,  we  developed  modified  versions  of  both  the  scaling  measures  and  the  previ¬ 
ously  published  entropy  measures.  These  extensions  include  the  generalization  of  the  entropy  as  well 
as  the  spatial  scaling  exponent,  S(h)  and  f(d)  respectively.  These  extensions  allow  us  to  not  only 
describe  the  average  seeding  measure  but  enables  us  to  assess  quantitatively  the  contribution  of 
different  subsets  to  the  overall  behavior.  This  approach  is  bas^  on  the  idea  that  the  animals  engage 
in  different  "types"  of  behavior  which  can  be  quantitatively  described  by  different  local  scaling 
exponents.  For  example,  a  rat  moving  along  the  wall  in  a  straight  path  may  also  engage  in  highly 
repetitive  head  movements  or  gnawing  stereotypies.  We  have  observed  precisely  this  type  of 
behavior  in  rats  treated  with  20.0  mg/kg  cocaine.  However,  the  averaging  proc^ure  yielding  the 
average  spatial  scaling  exponent  d  would  not  capture  the  distinctly  different  characteristics  of  these 
behavioral  "subsets".  Thus,  using  the  spectrum  of  spatial  scaling  exponents, we  can  now  assess 
these  behavioral  subsets  in  addition  to  the  average  geometrical  path  structure  [4].  The  singularity 
spectrum  of  local  entropies  ($$'’('  alpha  '  )  $)  computed  from  the  microcanonical  dynamical  par¬ 
tition  function  was  used  to  describe  the  sequential  dynamical  structure  of  rat  locomotor  paths  under 
the  influence  of  cocaine.  The  approach  us^  was  analogous  to  the  microcanonical  ensemble  of  sta¬ 
tistical  mechanics  and  permitted  the  sensitive  detection  of  "phase  transition"  like  behavior.  The 
results  suggest  that  the  behavioral  sequences  exhibited  by  the  rats  are  determined  by  two  processes, 
each  with  its  own  entropy  $  (  '  alpha  sub  (1)  ~  ,  "alpha  sub  (h)  ~  )  $.  These  processes  are  changed 
simultaneously  by  cocaine,  suggesting  that  they  coexist  as  local  repellers  with  overlapping  insets  and 
undergo  a  change  from  second  to  first  order  phase  transition  [8].  With  this  approach  we  may  be  able 
to  differentiate  between  different  drugs,  e.g.  MDMA,  high  doses  of  cocaine,  or  scopolamine,  gennat- 
ing  an  overall  increase  in  straight  path  patterns  but  acting  presumably  via  very  different  mechanisms. 
These  elaborations  have  already  1»1  to  even  more  precise  insights  into  mechanisms  underlying  the 
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behavioral  effects  of  psychoactive  drugs. 
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